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1. Introduction 

Random walks are a long studied problem not only in classical physics but also in many 
other branches of science [1]. Random walks help to understand complex systems, 
their dynamics and the connection between dynamics and the underlying topological 
structure. More recently, quantum analogues [2,3] of the classical random walks 
attracted considerable attention. 

Quantum walks have been introduced in the early nineties by Aharonov, Davidovich 
and Zagury [2]. Since that the topic attracted considerable interest. The continuing 
attraction of the topic can be traced back to at least two reasons. First of all, the 
quantum walk topic, of sufficient interest on its own right because there are 

fundamental differences compared to the classical random walk. Next, quantum walks 
offer quite a number of possible applications. One of the best known is the link between 
quantum walks and quantum search algorithms which are superior to their classical 
counterparts [4, 5] . Another goal is to find new, efficient quantum algorithms by studying 
various types of quantum walks. Let us point out that from a quantum mechanical point 
of view there is no randomness involved in the time evolution of a quantum walk. The 
evolution of the walker is determined completely by a unitary time evolution. This 
applies to the two most common forms: the discrete [2] and the continuous time [5] 
quantum walks. For a review on quantum walks see for instance [3]. 

Several aspects of quantum walks have been analyzed. First of all quite much 
attention was paid to the explanation of the asymptotics of the unusual walker 
probability distribution using various approaches [6-9]. Next, attention was paid to the 
explanation of the unusual probability distribution as a wave phenomenon [7, 13]. This 
question is of particular interest as there are several proposals how to realize quantum 
walks [10-17], in particular proposals using optical elements as the basic blocks for the 
quantum walks [12, 13]. Among the first ones was the optical implementation of a Galton 
board [14]. An ion trap proposal as well as an neutral atom implementation was put 
forward few years later [15-17]. The neutral atom proposal lead to a real experiment in 
the year 2003 [18]. 

Apart from the studies of elementary random walks its generalizations have been 
studied [19]. Generalizations of random walks to higher dimensions [20] have been put 
forward and differences to the simpler models pointed out. Next the effect of randomness 
in optical implementations on quantum walks was analyzed and its link to localization 
pointed out [21]. Among further generalizations also the behavior of more than one 
walker (particle) in networks realizing random walks has been studied [22,23]. The aim 
of the present paper is to add to this line of studies. We wish to study the evolution 
of two walkers performing a random walk. The evolution of each of the two walkers 
is subjected to the same rules. One of the interesting questions when two walkers are 
involved is to clarify how the probability of the walkers to meet changes with time (or 
number of steps taken in walk) . Because the single quantum walker behavior differs from 
its classical counterpart it has to be expected that the same will apply to the situation 
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when two walkers are involved. Interference, responsible for the unusual behavior of the 
single walker should play also a considerable role when two walkers will be involved. The 
possibility to change the input states (in particular the possibility to choose entangled 
initial walker states) adds another interesting point to the analysis. In the following, 
we will study the evolution of the meeting probability for two walkers. We point out 
the differences to the classical case and discuss the influence of the input state on this 
probability. 

The paper is organized as follows - first we make a brief review of the concept of 
the discrete time quantum walk on an infinite line and its properties. Based on this we 
generalize the two particle random walk for both distinguishable and indistinguishable 
walkers and define the meeting problem. This problem is analyzed in section 3. 
The asymptotic behavior of the meeting probability is derived. Further the effect of 
entanglement and indistinguishability of the walkers is examined. In the conclusions 
we summarize the obtained results and discuss possible development. Finally, in the 
appendix the properties of the meeting problem in the classical random walk are derived. 



2. Description of the walk 

We first briefiy summarize the description of the quantum random walk, for more details 
see e.g. [3]. We consider a coined random walk on an infinite line. The Hilbcrt space 
of the particle Ti. consists of the position space Tip with the basis {|m)|m G Z}, which 
is augmented by a coin space 

nc = {\L),\R)}. (1) 

The particle moves on the grid in discrete time steps in dependence on the coin state, 
the operator which induces a single displacement has the form 

S = ^\m-l){m\^ \L){L\ + ^ |m + l)(m| ® \R){R\. (2) 

m m 

A single step of the particle consists of a rotation of the coin state given by an arbitrary 
unitary matrix C and the conditional displacement S. The time evolution operator U 
describing one step of the random walk takes the form 

U^S{I^C). (3) 

If the initial state of the walker is |'?/'(0)), then after t steps its state will be given by the 
successive application of U on the initial state 

\m) = u'\m)- (4) 

The probability distribution generated by such a random walk is given by 

P{m,t) = \{m,L\i;{t))\' + \{m,R\m)\'- (5) 
In our paper we concentrate on a particular case when the rotation of the coin is given 
by the Hadamard transformation 
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since this is the most studied example of an unbiased random walk. For quantum walks 
with arbitrary unitary coin see e.g. [24]. 

We will describe the wave-function of the walker at time t by the set of two- 
component vectors of the form 

^(m,t) (7) 

where iljL{m,t) (^^(m, t)) is the probability amplitude that the walker is at time t on 
the site m with the coin state \L) {\R)). The wave-function thus has the form 

li^it)) = E (Mm, t) \m, L) + V'ij(m, t) |m, R)) . (8) 

m 

Throughout the text we will use symbols \'ip^^^\t)), ip^j^\m, t) for the vectors and 
the probability amplitudes, under the assumption that the initial state of the walker 
was 

|V^^(«)(0)) = |0, L(i?)), (9) 

similarly P^(^) (m, t) will be the corresponding probabilities. 

Prom the time evolution of the wave-function (4) follows the dynamics of the two- 
component vectors ■0(m, t) 

V'(m,i + l) = -^(^° + J J j^(m + l,i).(10) 

Thus the description of the time evolution of the walker reduces to a set of difference 
equations. Nayak and Vishwanath in [6] have found the analytical solution of (10) and 
derived the asymptotic form of the probability distribution. Before we proceed with the 
generalization of the quantum walk for two particles we will summarize the properties 
of the single walker probability distribution derived in [6], which we will use later for 
the estimation of the meeting probability. According to [6] the probability distribution 
of one walker is almost uniform in the interval (s — t/y^, s + ^/v^), where s is the initial 
position of the walker, with the value ^ 1/t. Around the points s ± t/\/2 are peaks 
of width of the order of t^/^ and the probability is approximately t~'^l^. Outside this 
region the probability decays faster then any inverse polynomial in t and therefore we 
will neglect this part. We will also use the slow oscillating part of the walker probability 
distribution derived in [6] which has the form 

Psli^{^.t) = —, 1 I (11) 

^^(l±f)^/l^ 

for the case that the initial coin state was L{R). 

The main difference between classical and quantum random walk on a line is the 
shape of the probability distribution. Due to the interference effect the quantum walker 
spreads quadratically faster than the classical one and its variance goes like cr ~ in 
contrast with the result for the classical case a ~ \/t. 

To visualize the properties of the quantum walk we plot in figure 1 the probability 
distribution and the slow-varying estimation. The initial conditions of the coin were 
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Figure 1. Probability distribution of a quantum random walker after 100 steps and 
the slow-varying estimation. 



chosen to be +i\R)), which leads to an unbiased distribution. The slow-varying 

estimation for the symmetric probability distribution is given by 

^ (12) 



^slow ^ 



TTt f 1 - f 



1 2x2 



Using this description of the random walk we now proceed to the analysis of the random 
walk for two non-interacting distinguishable and indistinguishable particles. 



2.1. Distinguishable walkers 

The Hilbert space of the two walkers is given by a tensor product of the single walker 
spaces, i.e. 

'H = {np®nc)i®{np®'Hc)2. (is) 

Each walker has its own coin which determines his movement on the line. Since we 
assume that there is no interaction between the two walkers they evolve independently 
and the time evolution of the whole system is given by a tensor product of the single 
walker time evolution operators (3). We describe the state of the system by vectors 

/ ipLL{m,n,t) \ 



'ijj(m, n, t) — 



(14) 



where e.g. the component iljfiL{m,n,t) is the amplitude of the state where the first 
walker is on m with the internal state \R) and the second walker is on n with the 
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internal state \L). The state of the two walkers at time t is then given by 

|V'(^)) = Z] Aj{m,n,t)\m,i)i\n,j)2. (15) 

m,nij= R^L 

The conditional probability that the first walker is on a site m at time t, provided that 
the second walker is at the same time at site n, is defined by 

P(m,n,t)= \{rn,i\{n,3\m)?- E \M'>^^n,t)\\ (16) 

i,j=L,R i,j=L,R 

Note that if we would consider a single random walker but on a two dimensional lattice, 
with two independent Hadamard coins for each spatial dimension, (16) will give the 
probability distribution generated by such a two dimensional walk. This shows the 
relation between a one dimensional walk with two walkers and a two dimensional walk. 
The reduced probabihties for the first and the second walker are given by 

Pi(m,t) =^P(m,n,t), P2{n,t) = J2 Pirn,n,t). (17) 

n m 

We can rewrite them with the help of the reduced density operators of the corresponding 
walkers 

Pi{t)^Trj^,m)){mi (18) 

in the form 

Pi{m,t)=J2{m,M{t)\m,j). (19) 
j 

The dynamics of the two walkers is determined by the single walker motion. Since 
we can always decompose the initial state of the two walkers into a linear combination of 
a tensor product of a single walker states and because the time evolution is also given by 
a tensor product of two unitary operators, the shape of the state will remain unchanged. 
Thus we can fully describe the time evolution of the two random walkers with the help of 
the single walker wave-functions. A similar relation holds for the probability distribution 
(16). Moreover, in the particular case when the two walkers are initially in a factorized 
state 

IV'(O)) = lY.^u{m,0)\m,i)A ® I ^^2,(n,0)|n,j)2) , (20) 

\m,i j \nj / 

which translates into iljij{m,n,0) = 'i/'ij(m, 0)'?/'2j(n, 0), the probability distribution 
remains a product of single walker probability distributions 

P{m,n,t) = {\^iL{m,t)\^ + \il:,R{rn,t)\^){\iP2L{n,t)\^ + \il:2R{n,t)\^) 

^ P,{m,t)P2{n,t). (21) 

However, when the initial state of the two walkers is entangled 

IV'(O)) = e| (EKKO)|m,z)i) ® (E^2,(^>0)Kj)2) I, (22) 

a I \m,i I \n,j I I 
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the probability distribution cannot be expressed in terms of single walker distributions 
but probability amplitudes 

2 



Yi^Um,t)^P^,{n,t) 



(23) 



P(m, n, t) = 

i,j=L,R 

As we have shown above, when the two walkers are initially in a factorized state, 
the probability distribution (21) is a product of the reduced distributions and thus the 
positions of the walkers are independent. On the other hand, when they are entangled, 
the probability distribution differs from a product of the reduced probabilities and the 
positions of the walkers are correlated. However, notice that the correlations are present 
also in the classical random walk with two walkers, if we consider initial conditions of 
the following form 

P(m, n, 0) = ^ (m, 0)P^{n, 0). (24) 

a 

The difference between (23) and (24) is that in the quantum case we have probability 
amplitudes not probabilities. The effect of the quantum mechanical dynamics is the 
interference in (23). 

Let us now define the meeting problem. We ask for the probability that the two 
walkers will be detected at the position m after t steps. This probability is given by the 
norm of the vector ip{m, m, t) 

MD{m,t)= Y \Aj{m,m,t)\^ = P{m,m,t). (25) 

i,j=L,R 

As we have seen above for a particular case when the two walkers are initially in a 
factorized state of the form (20) this can be further simplified to the multiple of the 
probabilities that the individual walkers will reach the site. However, this not possible in 
the situation when the walkers are initially entangled (22). The entanglement introduced 
in the initial state of the walkers leads to the correlations between the walkers position 
and thus the meeting probability is no longer a product of the single walker probabilities. 

2.2. Indistinguishable walkers 

We now analyze the situation when the two walkers are indistinguishable. Because 
we work with indistinguishable particles we use the Fock space and creation operators, 
we use symbols a^(^^-^ for bosons and for fermions, e.g. aj^j) creates one bosonic 
particle at position m with the internal state The time evolution is now given by 
the transformation of the creation operators, e.g. for bosons a single displacement is 
described by the following transformation of the creation operators 

«(m,L) ' [^^\m-\,L) + «(m+l,i?)) 

^{m.,R) ^ ('^{m-l,L) ~ '^(m+l,R)) ' (26) 

similarly for the fermionic walkers. The dynamics is defined on a one-particle level. 
We will describe the state of the system by the same vectors ip{m,n,t) as for the 
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distinguishable walkers. The state of the two bosonic and fermionic walkers analogous 
to (15) for the distinguishable is given by 

\^B{t))^Y. E ^(^iiK^'^) +V'ji(^,^,^))a(™,i)a{nj>«c), 

m,n iJ=L,R 

\Mt))^Yl Yl ^(V'ij(^,^,^) - V'ji(^,^,^))&(„,,i)&(„j)|'t^ac), (27) 

in,ni.j=L,R 

where \vac) is the vacuum state. In the case of two bosonic walkers initially in the same 
state we have to include an additional factor of 1/ ^/2, to ensure proper normalization. 
The conditional probability distribution is given by 

PB,F(m,n,t)= l(l(m,i)l(n,j)|V'B,F(t))r 
i,j=L,R 

= H \iJij{m,n,t) ±tljji{n,m,t)f, (28) 

i,j=L,R 

for m ^ n, and for m = n 

PB{m,m,t) = \{2^rn,L)\Mt))\''+\{'^{rn,R)\Mt))\'' + 
+ \{l(m,L)lim,R)\Mt))\'' 

= 2\jjjLL(m,m,t)\^ + 2\^jJRR{m,m,t)\'^ + 

+ li^LRim, m, t) + ipRL{m, m, t) |^ 
= MB{m,t), 

PF{m,m,t) = \{l(m,L)lim,R)\'4'Fit))\'^ 

= \'ipLR{m,m,t) -ipRL{m,m,t)\'^ 

= MF{m,t). (29) 

The diagonal terms of the probability distribution (29) define the meeting probability 
we wish to analyze. 

Let us now specify the meeting probability for the case when the probability 
amplitudes can be written in a factorized form ipij{m,n,t) — ipii{m,t)ip2j{'n,t) , which 
for the distinguishable walkers corresponds to the situation when they are initially not 
correlated. In this case the meeting probabilities are given by 

MB{m, t) = 2\ijjiL{m, t)i^2L{m, i)|' + 2|V'iR(m, i)V'2ii(m, t)\' 

m^ihim, t)ijj2R{m, t) + tpiR{m, t)ijj2L{m, t)\^, (30) 

for bosons and 

Mpim^t) = \ipiL{'m,t)i)2R{m,t) - Vii?(m, t)V^2L("^, ^)P, (31) 

for fermions. We see that they differ from the formulas for the distinguishable 
walkers, except a particular case when the two bosons start in the same state, i. e. 
-01(771, 0) = ip2{rn, 0) = ip{m, 0) for all integers m. For this initial state we get 

MB{m,t) = \Mm,t)t + \i}R{m,t)\' + 2\i}L{ni,t)ijR{m,t)\^ 
= {\ijL{m,t)\^ + \^PR{m,t)\^f 

= P\m,t), (32) 
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which is the same as for the case of distinguishable walkers starting at the same point 
with the same internal state. 

We conclude this section by emphasizing that except the case of two bosons with 
the same initial state the probability of meeting differs from the case of distinguishable 
particles and we have to use the probability amplitudes, whereas for the distinguishable 
walkers we can reduce it to one-particle probabilities. 



3. Analysis of the meeting problem 

Let us now compare the meeting problem in the classical and quantum case. We will 
study the two following two probabilities: the total probability of meeting after t step 
have been performed defined by 

M{t) = ^M{m,t), (33) 

m 

and the overall probability of meeting during some period of steps T defined as 
T 

M(r) = 1 - n (1 - ^(0) ■ (34) 

t=i 

The total probability of meeting M{t) is the probability that the two walkers meet at 
time t anywhere on the lattice, the overall probability of meeting M{T) is the probability 
that they meet at least once anywhere on the lattice during the first T steps. 



3.1. Distinguishable walkers 



We first concentrate on the influence of the initial state on the meeting probability for 
the distinguishable walkers. We consider three situation, the walker will start localized 
with some initial distance 2d (for odd initial distance they can never meet, without loss 
of generahty we assume that the first starts at position zero and the second at position 
2d), with the coin states: first, right for the first walker and left for the second 

V'i?L(0,2d,0) = l, (35) 

second, the symmetric initial conditions l/\/2(|L) + i\R)) for both 

/ 1 \ 



V'(0,2d,0) = 



(36) 



V -1 / 



and third, left for the first walker and right for the second 

^iR(0,2d,0) = l. (37) 

In the first case the probability distributions of the walkers are biased to the right 
for the first walker, respectively to the left for the second, and thus the walkers are 
moving towards each other. In the second case the walkers mean positions remain 
unchanged, as for this initial condition the single walker probability distribution is 
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Figure 2. Time evolution of the total probability of meeting for the three types of 
initial states and the classical random walk with two walkers. The initial distance is 
set to 10 and 20 points. 



symmetric and unbiased. In the last case the walkers are moving away from each other 
as their probability distributions are biased to the left for the first one and to the right 
for the second. 

Let us now specify the meeting probabilities. Recalling the expressions for the 
probability distributions and P^, we can write the total meeting probabilities with 
the help of the relations (21), (25) and (33) as 

MRL(t, d)^Yl P^'i^^ t)P''{m - 2d, t) 

m 

MLR{t, d)=Y. P'^im, t)P''{m - 2d, t). (38) 

m 

Wc sec that the meeting probability is fully determined by the single walker probability 
distribution. 

The figure 2 shows the time evolution of the total probability of meeting for the 
three studied situations and compares it with the classical case. The initial distance is 
set to 10 and 20 points. The plot clearly shows the difference between the quantum and 
the classical case. 

We derive some of the properties of the classical meeting problem in the Appendix 
A. The main results are the following. The meeting probabihty can be estimated by 

1 d"^ 

Mci{t,d)^^exp{—-). (39) 
y/TTt t 

This function has a maximum for t — 2d^, the peak value is given approximately by 

M,i{2d',d)^-jL^. (40) 

In contrast to the classical walk, in the quantum case the meeting probability is 
oscillatory. The oscillations come from the oscillations of the single walker probability 
distribution. After some rapid oscillations in the beginning we get a periodic function 
with the characteristic period of about six steps, independent of the initial state. In 
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the quantum case the maximum of the meeting probabihty is reached sooner than in 
the classical case - the number of steps needed to hit the maximum is linear in the 
initial distance d. This can be understood from the shape of the walkers probability 
distribution. The maximum of the meeting probability is obtained when the peaks of 
the probability distribution of the first and second walker overlap. If the initial distance 
between the two walkers is 2d then the peaks will overlap approximately after y/2d steps. 
The value of the maximum depends on the choice of the initial state. 

Let us derive analytical formulas for the meeting probabilities in the quantum case. 
For t > y/2d we consider the slowly varying part of the walker probability distribution 
(11), (12) and estimate the sums in (38) by the integrals 



MRLit,d) 



Ms{t,d) 



^ dx 



-'t\J^^ (1 - f)(l + ^)0^;i-2(^ 
t 

dx 



V2 

t 

V2 



MLR{t,d)^-^ [ "^"^ , (41) 

which can be evaluated in terms of elliptic integrals. Notice that the integrals diverge 
for d = 0, i.e. for the case when the two walkers start at the same point. We will discuss 
this particular case later, for now we will suppose that d > 0. The formulas (41) can 
expressed in the form 

MRL{t, d) ^ F+{2{t -d){t-{4- 2V2)d)K{a) + 

+ V2{{t - (4 + 2V2)d){t - (4 - 2V2)d)U{b+\a) - t'^U{c+\a))} 

Ms{t, d) f« ^^^+^~ {16d(t^ - d^){t + (4 + 2V2)d){t - (4 - 2V2)d)K{a) + 

+ V2(t + (4 + 2V2)d) (t - (4 + 2V2)d) (t+(A- 2V2)d) x 
X (i - (4 - 2^/2)d)((^ + d)U{b+\a) + (t - d)U(b_\a)) - 
- V2t'^{{t + d){t + (4 + 2V2)d){t + (4 - 2^/2)d)Il{c+\a) + 
+ (i - d){t - (4 + 2\/2)d){t - (4 - 2^/2)rf)^(c_|a))} 



MLR{t, d) ^ F_{2(t + rf)(t + (4 + 2V2)d)K{a) - 

- V2{{t + (4 + 2\/2)d){t + (4 - 2V2)(i)n(6_|a) - 
-t^n(c_|a))} 



(42) 
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where K(a) is the complete elhptic integral of the first kind and n(a;|a), n(|/|a) are the 
complete elliptic integrals of the third kind (see e.g. [25], chapter 17), a, b±, c± and F± 
are given by 

2t 



Ti^d{t T d){t{2 + \/2) T Ad){t{2 - \/2) ^ Ad)) 



[I ± ^/i){t - ^d) 
d{V2T'2) 



^ {t{V2T'2)+Ad){t-V2d) 

y/2d{t{^/2±2)-Ad) ' ^ ' 

From the relations (42) we can estimate the peak value of the meeting probability 

2-3^2 



MRL{V2d, d) 
Ms{V2d, d) 



7r(i(18 - 13^2) 
2 

ird 



The peak value shows a 1/d dependence on the initial distance between the two walker 
in all three studied cases, similar to the classical situation (40). 

Let us now analyze the asymptotic behavior of the meeting probability. In the 
Appendix A we show that in the classical case the meeting probability can be estimated 

by 

M^,(t,d)^^{l-^). (45) 
\/nt t 

In the quantum case we begin with the observation that the coefficients at the highest 
power of t with the elliptic integrals of the third kind arc the same but with the opposite 
sign for n(6|a) and n(c|a). Moreover, b± and c± go like —t as t approaches infinity, and 
thus all of the 11 functions have the same asymptotic behavior. Due to the opposite sign 
for n(6|a) and n(c|a) the leading order terms cancel and the contribution from this part 
to the meeting probability is of higher order oi 1/t compared to the contribution from 
the complete elliptic integral of the first kind K(a). The asymptotic of the function 
K{a) is given by 

K{a) ^ \.JU__ (46) 

t 

Inserting this approximation into (42) we find that the leading order term for the meeting 
probability in all three studied situations is given by (up to the constant factor a) 



ln( 



Moit, d) ^ a^^^. (47) 
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Figure 3. Left: Comparison of the meeting probability with the elhptic function 
estimation and the leading order term. The initial distance is set to 20 points. Right: 
Comparison of the meeting probability for the symmetric initial condition and zero 
initial distance with the estimation. 



Therefore, the meeting probabihty decays faster in the quantum case and goes hke 
\n{t/d)/t compared to the classical case (45). However, the decay is not quadratically 
faster, as one could expect from the fact that the single walker probability distribution 
spreads quadratically faster in the quantum walk. The exponential peaks in the 
probability distribution of the quantum random walker slow down the decay. 

The above derived results holds for d > 0, i.e. the initial distance has to be non- 
zero. As we have mentioned before, ior d — the integrals (41) diverge, and therefore 
we cannot use this approach for the estimation of the meeting probability. There does 
not seem to be an easy analytic approach to the problem. However, from the numerical 
results, the estimation 

Moit) ^ (48) 

fits the data best (6 being a constant prefactor). 

For illustration we plot in figure 3 the meeting probability and the estimations on 
a long-time scale. In the first plot is the case Mrl with the initial distance 20 points, 
on the second plot we have Ms and the initial distance is zero. 

Let us now focus on the overall meeting probability defined by (34). In figure 4 
we plot the overall probability that the two walkers will meet during the first T — 100 
steps. On the first plot we present the difference between the three studied quantum 
situations, whereas the second plot, where the meeting probability is on the log scale, 
uncovers the difference between the quantum and classical random walk. In the log scale 
plot we can see that the overall meeting probability decays slower in the quantum case 
then in the classical case, up to to the initial distance of ^/2T. This can be understood 
by the shape and the time evolution of a single walker probability distribution. After 
t steps the maximums of the probability distribution are around the point s ± \/2t, 
where s is the initial starting point of the random walker. For t — 100 steps the peaks 
are around the points s ± 140. So when the two walkers are initially more then 140 
points away, the peaks do not overlap, and the probability of meeting is given by just 
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Figure 4. The overall probability of meeting for two distinguishable quantum and 
classical walker during first 100 steps as a function of the initial distance. Same plot on 
the log scale. Only the values for even points are plotted since for odd initial distance 
the walkers never meet. 




Figure 5. Width of the overall meeting probability as a function of the upper bound 
T. 



the tails of the single walker distributions, which have almost classical behavior. Prom 
the first plot we see that the overall meeting probability is broader in the quantum case 
compared to the classical, which drops down very fast. The numerical results in figure 
5 show that the width of the overall meeting probability grows linearly with the upper 
bound T, the slope depends on the choice of the initial coin state. On the other hand 
in the classical case the width grows like y/T. 

Let us now analyze the overall meeting probability on a long time scale. In 
the Appendix A we show that in the classical case the overall meeting probability is 
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approximately given by 

M^(T, d)^l- exp |^-2^exp(-^) j exp (2dEric{-^)^ . (49) 

In the quantum case we consider In (1 — M£)(T, o?)) (Md stands for all three particular 
quantum cases) and estimate it with the help of (42) by 

T T 

In (1 - Md{T, d)) = ^In (1 - MdH, d)) ^ -Y.Mn{t, d) 

t=d t=d 
T 

^ - J MD{t,d)dt. (50) 

V2d 

Therefore we can estimate the overall meeting probabihty for T > \pld by 

Mi)(r, d) fti 1 - exp - y Mc(t, d)dt . (51) 

\ ^/2d J 

The meeting probability in the quantum case (42) involves elliptic integrals in a rather 
complicated form. However, we can estimate how fast the overall meeting probability 
converges to 1 for a fixed initial distance. This is determined by the rate at which the 
integral in (51) diverges. Before we proceed notice that 

ad 

j MD{t,d)dt^C{a), (52) 

i.e. the integral does not depend on the initial distance d. For large t we can estimate 
the meeting probability by (47) and thus we can divide the integral in (51) into 

/"^ 2 a/2 /"Inf^^l 

M^(t, d)dt ^ J Moit, d)dt + j y ' dt 

\Pid -Hd 

«C(«) + ^W(MI), (53) 

for appropriately large ol. Therefore the exponent in (51) goes like — In^(T) for large 
T. On the other hand, in the classical case from the estimation (49) we obtain that the 
asymptotic behavior of the exponent is given —\fT. Comparing these two we conclude 
that the overall meeting probability converges faster to one in the classical case. 

To conclude, the quadratic speed-up of the width shown in figure 5 follows from 
the the quadratically faster spreading of the quantum random walk. On the other hand 
for large times the meeting probability decays faster in the quantum walk, which leads 
to the slower convergence of the overall meeting probability. 

3.2. Effect of the entanglement for distinguishable walkers 



We will now consider the case when the two distinguishable walkers are initially 
entangled. According to (23) the meeting probability is no longer given by the product 
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of single walker probability distributions. However, it can be described using single 
walker probability amplitudes. We consider the initial state of the following form 

|V'(0)) = |0,2d)®|x), (54) 

where |x) is one of the Bell states 

\^^) = ^{\LR) ±\RL)) , 

\(t>^) = ^{\LL)±\RR)). (55) 

Recalling the probability amplitudes ip^^^\m,t) we can write the probability distri- 
butions of the two random walkers (16) in the form 

P^± (m, n, t) = 

I Y: |V.V,^)V'f(n-2d,^)±V'f(m,^)V'^(n-2d,^)|^ 

I Y: \^^{m,t)^|;f{n-2d,t)±^|;^(m,t)^|;f(n-2d,t)\'. (56) 

^ i,j=L,R 

The meeting probabilities are given by the sum of the diagonal terms in (56) 
M^±{t,d) = 

E \i^-{rri,t)ijf{m-2d,t)±^P^{rn,t)ij}{rn-2d,t)A, (57) 

^ m [i,j=L,R J 
M^±{t,d) = 

E \i'-{m,t)^Pf{m-2d,t)±ij^{m,t)i,f{m-2d,t)A. (58) 

^ m [i,j=L,R J 

The reduced density operators for both coins are maximally mixed for all four bell states 
(55). Prom this fact follows that the reduced density operators of the walkers are 

Pi(t) = ^(l^"(t))(V^"(t)| + |V^^(t))(V^«(t)|) 

P2{t) = I (li^mii^m + \i^fit)){i^m\) , (59) 

where \ip!i^^\t)) are analogous to \i/j^^^\t)) but shifted by 2d, i.e. 

H''^''\t)) = W{i;i^''\m,t)\m,L) + ^^^^''\m,t)\m,R)) 



m 



1 

2 



I E {^l''^ - t) \m, L) + (m - 2d, t) |m, R)) . (60) 



The reduced probabilities are therefore 

Pi(m,t) = ^(P^(m,t) + PV,^)) 

P2{m,t) = Pi{m-2d,t), (61) 
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Figure 6. Comparison of the meeting probability for the initially entangled coins and 
the symmetric case. The difference in the meeting probability. 



which are symmetric and unbiased. Notice that the product of the reduced probabihties 
(61) gives the probabihty distribution of a symmetric case studied in the previous 
section. Therefore to catch the interference effect in the meeting problem we compare 
the random walks with entangled coin states (55) with the symmetric case Ms- The 
figure 6 shows the meeting probabihties and the difference M^ — Ms, the initial distance 
between the two walker was chosen to be 10 points. 

We see that the effect of the entanglement could be both positive or negative. 
Notice that 

M^-{t, d) - Ms{t, d)^ - {M^+{t, d) - Ms{t, d)) 

M^-(t, d) - Ms(t, d)^ - (M^+(t, d) - Ms{t, d)) , (62) 

so the effect of \ip') is opposite to 10"*") and \4>~) is opposite to l'^'*'). The main difference 
is around the point t ~ V^d, i.e., the point where for the factorized states the maximum 
of the meeting probabihty is reached. The peak value is nearly doubled for M^- (note 
that for Ms the peak value is given by (44), which for d — 10 0.063), but 

significantly reduced for M^+ . On the long time scale, however, the meeting probability 
M^- decays faster than in other situations. According to the numerical simulations, the 
meeting probabilities for l-^"*") and \4>^) maintain the asymptotic behavior Int/t, but for 
it goes like 

M^-(t,d)^^. (63) 

V 

The initial entanglement between the walkers influences the height of the peaks giving 
the maximum meeting probability and affects also the meeting probability on the long 
time scale. 

3.3. Indistinguishable walkers 

Let us briefly comment on the effect of the indistinguishability of the walkers on the 
meeting probability. As an example, we consider the initial state of the walkers of the 
form |l(o,R)l(2(i,L))) i-G. one walker starts at the site zero with the right coin state and 
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Figure 7. Comparison of the meeting probability for bosons, fermions and 
distinguishable walkers. The difference in the meeting probability. 

one starts at 2d with left state. This corresponds to the case Mrl for the distinguishable 
walkers. The meeting probabilities are according to (30), (31) given by 

MB{t,d) = 5:{2|V^f (m,i)|>i(m - 2d,t)\' + 

m 

+ 2|V^f(m,t)nV'^(m-2d,i)|2 + 

+ l^f (m, t)ip^{m - 2d, t) + (64) 
+ ^|J^{m,t)^|J^{m - 2d,t\'^} 
Mj.{t, d)=J2 {li'Lim, t)iP^{m - 2d, t) - ^|(m, t)iPi{m - 2d, t)\^) . 

m 

In figure 7 we plot the meeting probabilities and the difference Mb,f — -^rl- From 
the figure we infer that the peak value is in this case only slightly changed. Significant 
differences appear on the long time scale. The meeting probability is greater for bosons 
and smaller for fermions compared to the case of distinguishable walkers. This behavior 
can be understood by examining the asymptotic properties of the expressions (64). 
Numerical evidence indicates that the meeting probability for bosons has the asymptotic 
behavior of the form hi{t)/t. For fermions the decay of the meeting probability is faster 
having the form 

MF{t,d)r^-. (65) 

The fermion exclusion principle simply works against an enhancement of the meeting 
probability. 

4. Conclusions 

We have defined and analyzed the problem of meeting in the quantum walk on an 
infinite line with two random walkers. For distinguishable walkers we have derived 
analytical formulas for the meeting probability. The asymptotic behavior following from 
these results shows that the meeting probability decays faster but not quadratically 
faster than in the classical random walk. This results in the slower convergency of 
the overall meeting probability. We have studied the infiuence of the entanglement 
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and the indistinguishability of the walkers on the meeting probabihty. The influence is 
particularly visible for fermions and in the case of distinguishable walkers for the case 
of initial entangled singlet state. 

Let us briefly comment on the correspondence between a one dimensional walk with 
two random walkers and a two dimensional walk. As two dimensional walks have been 
studied by many authors the possibility of common coin (i.e. a coin which is not a tensor 
product) for both walkers has arisen. In the context of one dimensional walk with two 
walkers this would mean some kind of interaction between the two walkers. It would be 
of interest to find coins which would attract the walkers and thus lead to an increase 
of the meeting probability, or repulsive coins with the opposite effect. Such interaction 
would be of infinite range as the walk will be driven by the same coin independent of the 
walkers distance. We can also think about local interactions where the walk is driven by 
a common coin for both walkers only when the distance between the walkers is smaller 
than some constant or they are at the same lattice point. 
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Appendix A. The meeting problem in the classical random walk 

Let us define the meeting problem on the classical level. We assume two particles which 
in each step of the process can perform randomly a step to the left or to the right on a 
one dimensional lattice labeled by integers. Initial distance between the two walkers is 
2d, because for odd initial distance the two walkers never meet, due to the transitional 
invariance we can assume that one walker starts in the vertex labeled by and the other 
one in the vertex 2d. We assume complete randomness, i.e. the probabilities for the 
step right or left are equal. We ask for the probability that the two particles meet again 
after t steps either at a certain position m or we might ask for the total probability to 
meet (the sum of probabilities at all of the possible positions) . Simple analysis reveals 
that the probability to meet at a certain position m equals to 

Mciit, m,d) = ^^ j ( ,+i-2d ) • (A-1) 

The total probability that the two particles are reunited after t steps reads 

Mci{t,d)= t ^il) (A.2) 

m=2d-t ^ \ 2 / \ 2 / 

which simplifies to 

This function has a maximum for t — 2d^ oi value 

(4d^V 

MJ2d^,d) = „ , ; „ , TTT. (A.4) 

''^ ' ' W\d{2d-l)y.{d{2d+l)y. ^ ' 

To obtain the asymptotic behavior of the meeting probability we approximate the one- 
walker probability distribution by a gaussian 

P,iix,t,d) = ^cxp(-^^^^), (A.5) 
Vnt 2t 



which leads to the following estimate on the meeting probability 



Y 1 d^ 

Ma{t,d)^ Pa{x,t,0)Pa{x,t,d)dx= ^eM-^)- (A.6) 
J \/Trt t 



-oo 



With the help of this estimation we can simplify the maximal probability of meeting 
into 

_ _ , .n 1 ... 
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For a fixed initial distance d we get the long-time approximation for t > d^ 

M^i{t,d)^^{l-^). (A.8) 

Finally we study the overall probability that the two walkers will meet at least once 
during first T steps, which is given by 

WiiT, 0?) = 1 - n (1 - Ma{k, d)). (A.9) 

k=d 



To estimate this function we take the logarithm of 1 — (T", d) and use the first order 
of the Taylor expansion to obtain 

T T 

/n(l - M^{T, d)) = J2 - Mciik, d))^-J2 Mdik, d). (A.IO) 

k=d k=d 

From this with the help of the relation (A. 3) we get the approximation of the overall 
meeting probability 

35;(r,d).i-axp(-|;^(/;j). (A.11) 

We estimate the sum in the exponential (A. 11) by an integral with the help of the 
formula (A. 6) 

S^l ( 2k\ J 1 d\ 

^exp(-^)-2dErfc(-^), (A.12) 
where Erfc is the complementary error function 

Eric{x) = 1 - Erf(x) = 1 - - / e'^'dt. (A. 13) 

TT J 


With the help of the estimation (A.12) we obtain the approximation of the overall 
meeting probability 

JU{T, d)^l- exp |^-2^exp(-^) j exp (2dErk{-^)^ . (A.14) 



